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polarization  from  bubbles  and  glass  spheres  in  water,”  Ocean 
Optics  VII  (Monterey,  June  1984). 
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VII  (Monterey,  June  1984). 

5.  P.  L*  Marston  and  D.  S.  Langley,  "Transmitted  wave  and  rainbow- 
enhanced  glories  of  dielectric  spheres,"  presented  at  the  1984 
Chemical  Research  and  Development  Center  Scientific  Conference  on 
Obscuration  and  Aerosol  Research  (Aberdeen  Proving  Ground,  MD, 

June  1984 )  . 

6.  P.  L.  Marston  and  J.  H.  Crichton,  "Radiation  torque  on  a  sphere 
illuminated  with  circularly  polarized  light,”  presented  at  the 
same  CRDC  Conference  as  the  preceding  item. 

The  following  papers  are  to  be  presented  at  the  108th  Meeting  of  the 
Acoustical  Society  of  America  (Minneapolis,  October  1984): 

7.  K.  L.  Williams  and  P.  L.  Marston,  "Axially- focused  (glory) 
scattering  due  to  surface  waves  on  tungsten  carbide  spheres." 

8.  P.  L.  Marston  and  D.  S.  Langley,  "Optical  scattering  properties  of 
bubbles  of  interest  in  acoustics  and  cavitation  research." 

9.  P.  L.  Marston  and  E.  H.  Trinh,  "Light  scattering  from  acoustically 
levitated  spheroidal  drops:  Discovery  of  a  new  hyperbolic  umbilic 
diffraction  catastrophe." 

C.  Technical  Reports  issued  under  this  Contract: 

1.  P.  L.  Marston,  Technical  Report  No.  3.  Focused  Acoustical  and 


Optical  Backscattering  from  Spheres  and  Properties  of  Qptoacoustic 


Sources,  Report  Number  AD-A133744  (Defense  Technical  Information 


Center,  Alexandria,  VA,  issued  September  1983). 

D.  Ph.D.  dissertations  supported  by  this  contract. 

1.  D.  S.  Langley,  "Light  Scattering  from  Bubbles  in  Liquids"  (Physics 
Department,  Washington  State  University,  July  1984). 

E.  Masters-Degree  Project  papers  (Note:  Washington  State  University 
allows  for  MS  degree  programs  in  physics  and  electrical  engineering 
which  have  no  requirement  of  a  formal  thesis.  The  students  in  these 
programs  are  required  to  submit  a  Project  Paper  which  is  reviewed  by 
faculty  members.  Normally  these  papers  are  not  distributed 
externally;  however,  they  sometimes  are  the  preliminary  description  of 
work  which  is  eventually  published.) 

1.  B.  L.  Brim,  "Theoretical  investigation  of  Brewster  angle  phenomena 
for  curved  surfaces"  (Dept,  of  Electrical  Engineering,  Washington 
State  University,  August  1983). 

II.  Preface  and  Explanation  of  System  of  References 

This  report  summarizes  progress  in  research  supported  by  the  contract 
titled:  "Propagation  and  Effects  of  Acoustical  and  Optical  Waves."  The 

emphasis  of  this  report  is  on  progress  subsequent  to  that  described  in  the 
report  by  P.  L.  Marston  titled: 

Technical  Report  No.  3.  Focused  Acoustical  and  Optical  Backscattering  from 
Spheres  and  Properties  of  Optoacoustic  Sources  (DTIC  No.  AD-A133744,  issued 
September,  1983). 

However,  for  continuity,  certain  research  discussed  in  that  report  will  be 


mentioned . 


The  following  reference  system  is  used  in  this  report.  References  to 

recent  external  communications  supported  by  this  contract  will  be  made  by 

giving  the  section  letter  and  number  of  the  list  given  in  Sec*  I  of  the 

present  report.  For  example,  the  first  item  listed  in  Sec*  I  is  referenced  as 
A1 

.  Reference  to  other  literature,  including  earlier  work  supported  by  this 

contract,  are  listed  in  Sec.  VII.  The  first  item  in  that  list  is  referenced 
1 

as 

Graduate  students  contributing  to  the  work  discussed  in  this  report 
include  Dean  S.  Langley,  Bruce  T.  Unger,  and  Kevin  L.  Williams.  The  research 
of  former  graduate  students  including:  Bradley  L.  Brim,  Timothy  J.  B.  Hanson, 
and  Jeffrey  L.  Johnson,  was  either  directly  used  and/or  will  merit  brief 
description  in  this  report.  Colleagues  at  institutions  other  than  Washington 
State  University  who  have  contributed  to  aspects  of  this  research  (without 
cost  to  the  contract)  include:  James  H.  Crichton  (Professor  of  Physics, 
Seattle  Pacific  University,  Seattle,  WA)  and  Eugene  H.  Trinh  (Staff  Scientist, 
Jet  Propulsion  Laboratory,  Pasadena,  CA) . 

III.  Acoustical  Scattering  Theory  and  Experiments 
A.  Motivation  for  this  Research 

The  motivations  for  this  study  of  the  scattering  from  elastic  objects 

in  water  are  as  follows:  (i)  It  is  an  example  of  the  development  and  testing 

of  new  models  for  the  scattering  from  elastic  objects  for  cases  in  which 

diffraction  provides  essential  corrections  to  elementary  ray  acoustics.  Such 

diffractive  corrections  occur  both  with  spherical  and  nonspherical  objects 

where  the  scattering  is  strongest  and  in  the  reflections  from  extended  curved 

surfaces.  (ii)  Our  models  consider  focused  backscattering  from  special 

spheres  which  should  have  giant  backscattering  cross  sections.^  Such  spheres 

2 

should  be  useful  as  calibration  targets  or  as  passive  navigational  beacons. 
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(Ill)  Our  models  consider  focused  Rayleigh  wave  contributions  to  the 
scattering  of  elastic  objects.  Such  contributions  are  particularly  indicative 
of  the  size  and/or  composition  of  the  scatterer*  (iv)  We  are  combining  our 
modeling  techniques  (which  use  geometrical  acoustics  and  diffraction  theory) 
with  other  developments  in  scattering  theory  including  "Resonance  Scattering 
Theory." 

One  salient  aspect  of  this  research  is  our  ability  to  model  and  measure 
scattering  amplitudes  not  just  in  backward  directions  but  also  in  near 
backward  directions*  It  is  shown  that  the  diffraction  lobes  in  near-backward 
directions  are  indicative  of  the  target’s  size. 

B*  Focused  backscattering  from  large  elastic  spheres: 

Transmitted  wave  acoustical  glories 

In  earlier  work,  we  have  shown  that  elementary  geometrical  acoustics 

predicts  an  infinite  backscattered  amplitude  due  to  certain  waves  transmitted 

through  a  sphere  but  that  diffraction  theory  may  be  used  to  remove  this 

unphysical  infinity.^  The  correctness  of  our  model  for  certain  classes  of 

glory  waves  was  experimentally  confirmed  using  the  apparatus  diagrammed  in 

A1 

Fig.  1.  The  first  confirmation  came  from  experiments  which  measured  the 
backscattering  due  to  unmixed-chord  glory  rays  such  as  those  diagrammed  in 
Fig.  2  for  a  sphere  of  fused  silica  in  water.  Figure  3  is  a  representative 
comparison  between  measurements  and  theoretical  predictions  when  ka  -  457 
where  a  is  the  sphere’s  radius,  k  =  2i t /X ,  and  X  is  the  acoustic  wavelength 
in  water.  The  spacing  of  the  diffraction  maxima  and  minima  in  near  backward 
directions  are  indicative  of  the  size  and  composition  of  the  scatterer. 

We  have  recently  extended  both  the  model  and  the  experiments  to  the 
much  more  difficult  problem  of  the  "mixed-mode"  or  "mixed-chord"  acoustical 
glory  from  an  elastic  sphere  in  water. ^  Our  previous  model^  was  complete 
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Fig.  1.  Simplified  diagram  of  the  scattering  experiments.  The  hydrophone 
be  scanned  along  a  line  transverse  to  the  symmetry  axis  defined  by  the  sou 
and  the  sphere.  The  diagram  is  not  drawn  to  scale. 
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Fig.  2.  Representative  unmixed-chord  glory  rays  for  a  fused-silica  sphere  in 
water..  Dashed  and  solid  rays  represent  shear  and  longitudinal  waves 
respectively. 
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Fig.  3.  The  peak-to-peak  voltage  for  the  (4,4,1)  echo  plotted  as  a  function 
of  x.  The  voltages  are  ratioed  with  the  on-axis  value  of  the  first  axial 
echo.  Measurements  were  obtained  when  the  water  temperature  was  24.9°  C.  The 
solid  curve  is  the  modeled  dependence  on  x  =  0.  The  peak  at  x  =  0  and  the 
side  lobes  are  direct  evidence  of  axial  focusing. 
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Table  I.  Predicted  and  Measured  On-Axis  Pressure  Ratios  for  the  (3,1,1)  and  (3,2,1)  Echoes  at  Two  Dif¬ 
ferent  Receiver-Sphere  Distances.  Also  shown  is  B  and  Y  for  each  chord  sequence  at  both 
echoes.  Tjie  water  temperature  was  20°C.  °  a 
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only  for  transmitted  waves  within  the  sphere  which  were  either  all 
longitudinal  (L)  or  all  shear  (S).  In  this  new  research  we  completed  the 
theory  for  the  case  of  mixed-mode  echoes  in  which  there  is  at  least  one  L-to- 
S  or  S-to-L  mode  conversion  due  to  internal  reflection*  The  chord  sequence 
of  one  of  the  mixed  modes  studied  is  shown  in  Fig.  4.  A  novel  calculation 
shows  that  the  focal  parameters  are  independent  of  permutations  of  the  mode 
sequence;  however,  the  contributions  to  the  scattering  depend  on  sequence. 

That  a  phase  factor  and  an  amplitude  factor  B  depend  on  the  sequence 

is  shown  by  the  theoretical  values  tabulated  in  Table  I.  The  full  expressions 
for  the  amplitude  and  phase  are  given  in  Ref.  A8.  For  the  purposes  of  the 
present  discussion  it  is  Important  to  note  that  a  correct  prediction  for  the 
mixed-mode  amplitude  is  only  possible  If  the  interference  between  the 
different  sequences  contributing  to  that  mode  are  properly  modeled. 

Data  which  confirms  our  model  of  the  mixed— mode  backscattering  is 
displayed  in  Table  I.  and  are  experimental  and  theoretical 

amplitudes,  respectively,  which  are  normalized  to  the  amplitude  of  the  first 
specular  reflection  from  the  sphere.  In  some  cases  the  agreement  between 
experiment  and  theory  is  seen  to  be  better  than  1%.  Data  acquisition  and 
analysis  was  made  possible  by  using  a  digital  signal  averager  and  a  digital 
tone  burst  generator  both  purchased  with  contract  funds. 

G.  Axially— focused  (glory)  backscattering  due  to  Rayleigh  waves 
generated  on  elastic  spheres. 

For  most  sonar  systems  and  targets,  the  ka  is  significantly  smaller 
than  the  ka  of  the  aforementioned  experiments  on  transmitted-wave  glories.  It 
may  be  argued  that  when  ka  <  100  certain  radiation  damped  "surface"  waves  will 
significantly  contribute  to  the  backscattering  from  elastic  spheres  or 
spherical  shells.  Which  type  of  surface  wave,  and  their  relative  strengths, 
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will  depend  on  the  physical  nature  of  the  scatterer.  It  is  of  general  interest 
both  to  demonstrate  (and  to  model)  the  axial— focusing  of  backscattering  due  to 
surface  waves.  We  have  measured  and  modeled  the  particular  case  of  scattering 
from  tungsten  carbide  spheres  in  water  in  the  range:  30  £  ka  _<  100.  The 
dominant  surface  wave  in  this  range  is  the  Rayleigh- like  wave.  The  angular 
dependence  of  its  contribution  to  the  scattering  was  modeled  and  the  model  was 

B7 

confirmed  with  measurements.  The  angular  dependence  is  indicative  of  the 
sphere's  size. 

The  experiments  were  carried  out  with  the  measurement  configuration 
diagrammed  in  Fig .  1 .  The  apparatus  was  changed  ( from  that  used  in  our 
previous  work  *  )  to  extend  its  bandwidth  and  reduce  its  operating  frequency 

by  changing  the  transducers  and  the  preamplifier.  The  tungsten  carbide 
spheres  had  diameters  of  1  inch  and  0.5  inches.  When  a  short  tone  burst  is 
incident  on  a  tungsten  carbide  sphere,  the  backscattering  response  shows  a 
decaying  pulse  train  which  can  be  attributed  to  a  specular  reflection  and  the 
repeated  circumnavigations  of  waves.  Representative  cases,  measured  with  our 
apparatus,  are  shown  in  Fig.  5.  The  specular  and  first-Rayleigh-wave 
contributions  to  the  backscattering  are  identified.  The  general  appearance  of 
these  curves  are  similar  to  a  measurement  of  the  transient  backscattering 
published  by  the  group  at  the  Naval  Research  Laboratory.^  In  the  NRL 
measurement,  however,  the  incident  sound  was  not  a  well  defined  burst  of  sine 
waves  (as  in  Fig.  5)  and  the  angular  dependence  of  the  scattering  (discussed 
below)  was  not  described. 

If  the  receiver  is  moved  away  from  the  backscattering  axis,  the 
amplitude  of  the  Rayleigh  wave  contribution  at  first  decreases,  and  then 
undergoes  an  oscillatory  dependence  on  the  transducer  displacement  x.  We  have 
developed  a  theory  for  the  angular  dependence  which  predicts  the  Rayleigh  wave 
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contribution  to  the  scattered  sound  has  a  pressure  amplitude  which  satisfies 
the  following  proportionality: 


fWleigh(ka*Y) 


"I, 


PRa7leigh<ka-Y-°)J0CSY) 


(1) 


where  Y  =  tan  (x/(z+a))  is  the  angle  relative  to  the  backscattering  axis, 

«  1  radian,  and  is  a  Bessel  function.  The  constant  3  is  predicted  to 

depend  on  the  frequency  of  the  sinusoidal  tone  burst  used  (or  equivalently  the 
ka  of  the  sphere).  The  general  reasoning  leading  to  Eq.  (1)  is  similar  to 
that  which  we  give  in  References  1,  Al,  and  A5  for  other  cases  of  axially- 
focused  scattering:  the  scattered  pressure  wave  in  the  water  appears  to 
emanate  from  a  virtual  ring- like  source  which  we  call  the  focal  circle. 

We  first  tested  the  proportionality  on  J^(3y)  by  measuring  the 
dependence  of  Rayleigh  amplitude  and  plotting  the  ratio 
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^Rayleigh 

1 PRayleigh  (ka,Y=0> 


(2) 


Representative  data  points  for  the  case  ka  =  59.2  are  shown  in  Fig.  6.  An 

empirical  value  for  3>  call  it  was  determined  by  adjusting  J  (3  Y)  to 

"  0  E 

best  fit  these  data.  This  procedure  is  also  illustrated  in  Fig.  6  which  shows 

that  these  data  are  of  the  expected  form.  Using  this  empirical  3  ,  Eq.  (1) 

E 

predicts  the  first  null (where  PRayleighCka»^1)  =  0)  to  be  at: 


2.4048 

Yt  -  — 7, -  radians  . 

i  Pt? 


(3) 


We  have  also  developed  a  theory  for  3  which  is  based  on  finding  the 
focal-circle  radius  from  a  phase  matching  condition.  The  result  is 


ka  c/c 

P  * 
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where  c  is  the  velocity  of  sound  in  water  and  is  the  phase  velocity  of 

the  Rayleigh  wave  on  a  tungsten-carbide  sphere  in  water  determined  at  the  ka 
of  the  incident  tone  burst.  The  procedure  for  calculating  c/c  is 

—  p 

complicated  and  will  only  be  briefly  summarized.  An  extension  of  the  theory 

of  resonance  scattering  (Ref.  4,  Eq.  46c)  gives  c/c  as 

P 


Re(X  ) 

n.l 


c 

c  n  +  (1/2)  5 


(5) 


at  discrete  values  ka  =  Re(X„  ,  ) .  Here  X^  ,  is  the  Rayleigh  complex  root 

I*  J  -L  II  y  J.  — — 1  "  ■■ 

of  the  following  equation 


Wi’  *  0  • 


(6) 


in  which  is  the  following  determinant 
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d22 

d23 
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(7) 

0 

d32 

d33 

and  the  d„(X)  in  this  determinant  are  the  complex  functions  listed  in  the 

Appendix  of  Ref.  5.  The  d„  are  functions  of  X  involving  spherical  Bessel 

and  Hankel  functions,  j  and  h  ,  their  derivatives,  and  the  elastic 

n  n 

parameters  of  the  sphere  and  of  the  surroundings  (water).  Equation  (6)  was 

solved  numerically  for  the  roots  X  for  the  partial-wave  index  n  having 

integer  values  from  4  to  30.  This  procedure  (carried  out  by  Kevin  Williams) 

was  facilitated  by  adapting  computer  programs  previously  implemented  by  B.  L. 

E 1 

Brim  during  previous  contract  related  research.  The  values  of  X  agreed 

n ,  1 

with  those  given  in  Table  IV  of  Ref.  5  for  the  specific  cases  n  =  4  to  7 . 

The  Cp/ c  values  thus  obtained  from  Eq.  (5)  for  discrete  ka  <  85  were  fitted 
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with  a  polynomial  to  give  c  /c  as  a  continuous  function  of  ka  over  the  ka- 

P 

range  of  interest. 

The  aforementioned  procedure,  together  with  Eqs.  (1)  and  (4),  gives  the 
following  theoretical  value  for  the  null  angle 


Y1 


2.4048  . V 
ka  1  c  ; 


radians  . 


(8) 


Figure  7  compares  the  theoretical  (y^)  ^  (solid  line)  with  the  experimental 
values  (circles)  from  Eq.  (3).  The  ordinate  is  the  ka  value  for  the  incident 
burst.  The  agreement  between  theory  and  experiment  is  sufficient  to  confirm 
the  principal  features  of  our  model  of  the  angular  width  of  the  focused 
backscattering . 

There  may  be  a  small  systematic  source  of  error  in  the  theoretical 

calculation  displayed  in  Fig.  7  because  the  elastic  properties  of  the  specific 

tungsten  carbide  samples  used  in  this  experiment  were  not  precisely  known. 

Calculations  of  c  / c  based  only  on  tabulated  elastic  properties  were 

adjusted  to  according  to  the  measured  group  velocity  c  for  these  spheres 

g 

according  to  a  relation  between  cg  and  c^  derived  by  Williams  for  the  ka 
region  of  interest.  The  resulting  theory  line  is  the  one  shown  in  Fig.  7.  If 
only  previously  tabulated  elastic  properties  are  used,  the  revised  theory  line 
has  a  slightly  larger  slope.  The  theory  lines  computed  with  and  without  the 
aforemention  correction  bracket  the  data. 

Though  the  experiments  we  described  here,  and  in  Sec.  IIIC,  use 
incident  tone  bursts  of  only  a  few  cycles  in  length,  the  concept  of  axial 
focusing  is  also  applicable  to  spheres  or  spherical  shells  insonified  by  long 
incident  tone  bursts.  For  example,  certain  large  resonance  features  in  the 
scattering  from  spherical  shells6  may  be  attributed  to  a  particular  partial 
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wave,  say  the  nth  partial  wave.  The  angular  dependence  of  the  amplitude  for 
the  nth  partial  wave  may  be  expanded  using  Szegofs  asymptotic  expansion  of  the 
Legendre  polynominal^ 


Pn(cosY)  ^  J0[(n  +  |)y3  , 


(9) 


when  n  »  1  and  y  «  1  radian.  This  is  of  the  same  form  we  have 
demonstrated  here  in  Fig.  3  and  6  for  selected  responses  to  short  tone  bursts. 

D.  Nearfield  effects  in  acoustic  scattering 

We  have  investigated  aspects  of  the  nearfield- to-farfield  transition 

for  the  backscattering  from  an  elastic  sphere.  The  motivation  for  this 

research  is  as  follows:  (i)  Following  a  presentation  (at  the  Acoustical 

7  A1 

Society  of  America  Meeting,  Orlando,  FL,  Nov.  1982)  of  our  research,  * 
questions  were  raised  by  a  member  of  the  audience  concerning  the  theoretical 
form  (Eq.  10,  below)  used  to  model  the  specular  (i.e.,  first  axial) 
contribution  to  the  scattering.  A  detailed  derivation  of  this  result  was 
subsequently  published  as  an  insignificant  part  of  Ref.  A1 .  However,  it  is 
not  too  surprising  that  our  model  of  the  specular  reflection  is  not  well  known 

since  it  differs  in  detail  with  that  used  in  the  NRL  study  of  calibration 

3  7 

spheres.  (ii)  There  has  been  recent  interest  in  the  criterion  for 

establishing  the  farfield  in  scale  model  scattering  experiments. 

With  reference  to  Fig.  1,  our  model  for  the  amplitude  of  the  first 

A1 

axial  (or  specular)  reflection  for  a  tone  burst  having  ka  »  1  is 


2(z  +  j  a) 


®WW 


-ctz 


(10) 


where:  a  is  the  sphere’s  radius;  |pinc(  is  the  pressure  magnitude  of  the 
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incident  tone  burst  measured  at  point  C'  in  Fig.  1;  a  is  the  attenuation 
coefficient  of  sound  in  water  for  the  central  frequency  of  the  tone  burst;  and 

*  (piCL  **  poCW)/(piCL  +  PoV*  Here  pi  and  po  (and  CL  and  V  arS 
respectively  the  densities  (and  longitudinal  sound  speeds)  of  the  elastic 

sphere  and  of  water.  The  feature  at  question  here  are  the  (z  + a) 
divergence  factor  which  shows  the  reflection  appears  to  emanate  from  a  point¬ 
like  source  located  half-way  between  C*  and  C  in  Fig.  1.  Also  important  is 
our  assertion  that  is  less  than  unity  in  Eq.  (10).  Our  derivation  of 

Eq.  (10)  assumes  not  only  that  ka  »  1  but  also  that  the  source  distance  z^ 

»  a  and  that  the  (small)  receiver  hydrophone  is  located  on  the  backward  axis 
at  x  =  0  in  Fig.  1. 

To  test  the  correctness  of  Eq.  (10),  measurements  were  made  of  the 
specular  amplitude  |p^|  at  several  distances  z  from  a  fused  silica  sphere 
having  a  =  51.6  mm.  The  z  ranged  from  6  cm  to  ~95  cm.  The  frequency  of 
the  incident  tone  burst  was  2  MHz.  Comparison  of  the  measured  |p0l/lpincl 
with  the  predictions  of  Eq.  (10)  clearly  demonstrate  that  (z  +  ^  a)  is  the 
correct  divergence  factor  in  Eq.  (10).  Indeed,  if  (z  +  ^a)  is  replaced  by 
(z  +  a)  ^  (which  is  the  form  often  assumed  in  farfield  models  of  the 
scattering),  the  modified  theory  quite  clearly  disagrees  with  the  data, 
especially  for  z  <  50  cm.  Note  that  (z  +  a)  is  the  distance  of  the 
hydrophone  from  the  sphere’s  center. 

1  9a 

In  related  developments,  we  have  previously  noted  ’  that  a  general 
farfield  condition  of 

z  »  ka2  (ID 

follows  from  the  requirement  that  the  farfield  is  the  region  where 
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Fraunhofer’s  approximation  to  diffraction  integrals  is  applicable*  We  have 
also  shown  Eq.  (11)  by  answering  the  question:  ’’When  are  farfield 
approximations  applicable  in  exact  partial-wave  theories  for  the  scattered 
field?”  Surprisingly  the  condition  (11)  does  not  appear  to  be  well  known 
since  the  principal  investigator  (P.  L.  Mars ton)  recently  refereed  a 
manuscript  concerned  with  nearfield  scattering  in  which  (11)  was  evidently  not 
known  to  the  author  even  though  his  calculations  supported  this  result!  At 
Marston’s  suggestion,  a  discussion  of  Eq.  (11)  included  and  the  manuscript  was 
published  as  Ref.  9b.  The  condition  Eq.  (11)  should  be  a  sufficient  condition 
in  acoustic  scattering  from  bounded  nonspherical  objects  provided  a  is  taken 
to  be  half  the  largest  width  of  the  object. 

E«  Backseat tering  from  a  large  torus:  a  simple  example  of 
axial  focusing 

The  scattering  from  a  large  torus  was  modeled  for  the  case  of  an 

incident  plane  wave  which  propagates  parallel  to  the  symmetry  axis  of  the 
A5 

torus.  The  configuration  modeled  is  shown  in  Fig.  8.  In  this  case  the 
specularly  reflected  wave  is  axially  focused;  it  appears  to  emanate  from  a 
focal  circle  (at  F^  in  Fig.  8).  For  a  large  rigid  torus  the  pressure  of  the 
reflected  wave  is  predicted  to  be^ 


p  #  b  „ 

=  ( — ~ — )  (irka)  J(/u0) 


i  (kr  +  Cp  -  —  -  o)t) 


(12) 


where:  =  kbsiny ;  a  and  b  are  as  shown  in  Fig.  8;  the  polar  coordinates 

(relative  to  CT)  of  the  observation  point  are  r  and  Y;  <f  ■  -ka(l  -  cosY)/2; 

an<^  Pine  amplitude  at  Cf  of  the  incident  plane  wave. 

This  model  was  developed  so  as  to  identify  an  example  of  focused 

backscattering  which  would  be  free  of  some  of  the  intricacies  of  the  modeled 

1  A1 

scattering  from  spheres.  9 


Observations  of  scattering  from  a  torus  would  be 


VIRTUAL 
WAVE 
FRONT  * 


Fig.  8.  Cross-sectional  view  of  the  rigid  torus  for  which  the  backscattering 
was  modeled.  The  figure  may  be  rotated  about  the  CC'  axis.  The  reflected 
wave  is  axially  focused  along  the  left  extension  of  the  CC'  axis  when  the 
incident  wave  propagates  in  the  -z  direction.  Scattering  due  to  the  dotted 
ray  is  omitted. 
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A5 

particularly  useful  in  testing  a  model  of  pulse  distortion  which  results 
from  focusing.  There  are,  however,  other  potential  applications  which 
include:  (i)  the  use  of  a  torus  as  either  a  calibration  target  or  as  a 

navigational  target  because  of  the  simple  form  manifest  in  Eq.  (12);  (ii)  the 
use  of  a  torus  as  a  special  case  in  which  to  test  computational  algorithms  for 
the  scattering  from  objects  having  large  aspect  ratios.  (There  have  been 
recent  advances ^  in  the  development  of  algorithms  applicable  to  prolate 
scatters  having  large  aspect  ratios  and  similar  advancements  for  the  case  of 
oblate  scatters  are  to  be  anticipated.)  Unless  the  conditions  ka  »  1,  kb  » 
1,  and  b  »  a  are  all  met,  however,  in  addition  to  (12)  there  will  be  other 
significant  contributions  to  the  backscattering  including  those  due  to  the 
dotted  ray  in  Fig.  8  and  to  surface  waves  which  loop  around  either  side  of  the 
torus . 


F.  Theory  of  pulse  distortion  in  focused  scattering:  half-order 
derivative  of  a  sine-wave  burst 

An  important  property  of  focused  scattering,  such  as  that  discussed  in 

Sec.  IIIB  and  IIIE  of  this  report,  is  its  large  magnitude.  Then  it  is 

reasonable  to  inquire  whether  certain  features  of  a  scattered  transient  are 

indicative  of  focused  scattering  If  such  features  were  unmistakable,  it  would 

be  possible  to  discriminate  between  sonar  targets  which  were  large  in  size 

(and  having  no  focused  backscattering)  and  targets  which  were  small  (whose 

signals  were  large  because  of  focusing  of  the  backscattering) .  Our 

calculation  of  the  half-order  derivative  of  a  sine-wave  burst^3  is  a  step 

toward  the  resolution  of  questions  of  this  type.  In  our  model  of  focused 

A1 

high-frequency  scattering  from  spheres,  it  was  found  that  certain  of  the 
axially  focused  waves  had  pressures  proportional  to 
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S(X)  [  (t-T)TT ] ~  1/2  dT  , 


(13a) 


where  s(t)  =  ds/dt  at  t  *  t  and  s(t)  is  proportional  to  the  pressure  of 
the  incident  wave  p  •  In  t^ie  Present  calctxlation  the  incident  wave  was 
modeled  as  a  burst  of  sine  waves  N  cycles  in  length  so  that  s(t)  could  be 
written  as  the  product 

s  (t)  =  H(co t)H(2TTN  -  o)t)sinojt  ,  (13b) 


where  the  Heaviside  unit-step  function 


H(y  <  0)  =  0,  H(y  >  0)  =  1,  and  H(y  «  0)  »  1/2. 

1  /2 

It  was  found  by  direct  evaluation  of  Eq.  (13)  that  (d/dt)  may 

be  expressed  using  a  combination  of  phase  advanced  sine  waves  and  the 

A5 

auxiliary  Fresnel- integral  function.  The  leading  edge  of  a  focused  burst  is 
steepened  and  the  peak  amplitude  of  the  first  half-cycle  is  suppressed.  The 
trailing  edge  of  the  focused  burst  is  predicted  to  manifest  a  tail.  Before 
recent  improvements  in  our  experimental  apparatus.  Fig.  1,  the  bandwidth  of 
our  system  was  judged  to  be  insufficient  to  experimentally  test  these 
predictions  in  detail.  Results  of  this  calculation  for  the  central  peak-to- 
peak  amplitude  of  a  focused  burst  having  N  >  4  justify  the  experimental 
procedures  used  in  the  confirmation  of  theory  in  Refs.  A1  and  A8. 

In  addition  to  aforementioned  motivations  for  the  evaluation  of  Eq. 
(13),  the  following  should  be  noted.  It  would  be  possible  to  compare  these 
results  with  those  of  an  "exact"  calculation  of  the  scattering  of  a  burst 
based  on  the  evaluation  of  the  partial-wave  series  at  many  discrete 
frequencies  and  the  use  of  a  fas t-Fourier- transform  algorithm.  Such  an  exact 
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computation  is  sufficiently  intricate  that  the  comparison  of  results  with 
simple  models  (for  special  cases)  is  justified. 

IV.  Scattering  of  Light  from  Bubbles  in  Liquids  and  Other  Problems  in 

Electromagnetic  Scattering  of  Value  to  Acoustics  and  Naval  Technology 
A.  Motivations  and  Overview 

The  majority  of  the  contract  research  into  electromagnetic  scattering 
is  concerned  with  the  optical  properties  of  bubbles  in  liquids.  In  this 
section  some  motivations  for  understanding  the  scattering  properties  of 
bubbles  will  be  noted. 

Gas  bubbles  of  natural,  origin  inhabit  the  upper  layers  of  the  ocean;^ 

however,  the  causes,  lifetimes,  and  size  spectra  of  these  natural  bubbles  are 

12 

not  well  understood.  Medwinfs  acoustic  measurements  suggest  that  the 
bubbles  may  be  present  in  large  numbers  even  at  depths  of  36  meters  where  he 
estimated  (Fig.  6  of  Ref.  11)  there  to  be,  for  example,  -  8000  bubbles/m  in 

the  radius  interval  65  to  66  ym.  Bubbles  also  affect  the  colors  of  the  sea  in 

13  1  a 

regions  of  breaking  waves,  and  the  optical  properties  of  sea  ice. 

Bubbles  are  important  to  various  naval  technologies.  For  example, 

natural  bubbles  may  result  in  fluctuations  of  the  acoustic  properties  of  the 

near-surface  ocean  (Ref.  11  and  papers  cited  therein).  Man  made  bubbles  may 

15  i  & 

be  injected  in  the  ocean  to  regionally  reduce  the  transmission  of  sound  ’ 

(the  bubble  screen"  technique.)  Microbubbles  also  govern  the  inception  of 

17 

acoustic  and  hydrodynamic  cavitation.  Microbubbles  alter  the  effective 
nonlinearity  parameter  (B/A)  of  a  liquid  and  also  the  generation  of  acoustic 
2nd  harmonics  and  parametric  signals  (having  the  sum  or  difference  frequency 
of  the  pump  waves) .  Hydrodynamic  flow  drag  may  be  reduced  by  the  injection  of 
microbubbles  (through  a  porous  section  of  the  wall)  into  the  liquid  boundary 
layers.^  Bubbles  are  also  useful  as  tracers  in  hydrodynamics  research 3® 


Bubbles  may  also  be  significant  in  optical  communications  and/or  detection 

technologies  involving  submerged  platforms. 

In  spite  of  the  aforementioned  significance  of  bubbles  (and  other 

applications  noted  below)  no  systematic  research  into  the  optical  scattering 

properties  of  bubbles  had  been  carried  out  prior  to  the  present  ONR  Contract 

research  which  began  in  1980.  Critiques  of  the  research  prior  to  about  1980 

are  given  in  Ref.  A4,  A9,  A10,  19,  and  20  and  we  need  only  summarize  them. 

21 

The  principal  modeling  approaches  had  been  based  on  geometrical  optics  and 

significant  effects  of  interference  and  diffraction  were  not  included.  Such 

models  might  appear  to  be  reasonable  since  bubbles  are  usually  much  larger 

than  the  wavelength  of  light  in  water;  we  find,  however,  that  geometrical 

methods  by  themselves  are  deficient  even  for  large  spherical  bubbles. 

Deficiencies  are  evident  in  the  various  transition  regions  of  the  scattered 

light,  particularly  in  the  near  forward  and  backward  directions  9  and  in 

19  A4 

the  critical-angle  region  *  where  the  transition  from  partial  to  total 

reflection  at  the  bubble1 s  surface  occurs*  In  the  light  of  our  Mie 
19 

calculations  (which  are  exact  for  spheres) ,  physical-optics 
23  A4  A10 

approximations,  and  (now)  data  *  the  early  (1972)  Mie  calculations  of  A. 
24 

Keller  appear  to  be  erroneous*  (We  mention  this  because  Keller ! s  research 
is  still  used  to  justify  certain  optical  technologies  for  counting  bubbles  in 
cavitation  research. ) 

In  the  present  research  we  continue  to  model  bubbles  as  spherical  voids 
surrounded  by  a  homogeneous  dielectric  media  (usually  water).  Some  recent 
data  support  this  assumption  for  small  bubbles  in  clean  or  moderately  lean 
water.  Even  with  the  simplification  of  an  assumed  spherical  shape,  the 
backscattering  properties  of  bubble- like  objects  were  not  considered  by 
previous  research  into  short-wavelength  scattering  (including,  e.g.  well  known 
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previous  research  into  short-wavelength  scattering  (including,  e.g.  well  known 

25  26  27 

papers  by  Rubinow,  Nussenzveig,  and  Inada  and  Plonus^ ) . 

There  is  already  evidence  that  the  present  line  of  research  has 

affected  existing  or  future  technologies  relating  to  the  measurement  of  bubble 

12  28 

size  spectra  in  the  ocean,  *  the  measurement  of  bubble  sizes  in  acoustics 

29 

cavitation  research,  and  the  modeling  of  bubbles  in  glass  (as  applied  to 

30 

high-power  laser  glass  amplifiers  ).  It  is  likely  to  affect  research  on 
bubble  nucleation  in  polymer  liquids  (as  applied  to  plastics  technology11). 

The  emphasis  of  this  research  has  been  on  the  understanding  of 
fundamental  scattering  properties  of  bubbles  as  noted  in  Sec.  IVB-IVD  of  this 
report.  These  properties  are  relevant  to  various  bubble  detection  and  sizing 
technologies.  We  also  model  a  novel  application  in  acoustics  in  Sec.  IVE: 
the  potential  use  of  light  scattering  to  directly  measure  the  response  of  a 
bubble  to  an  acoustic  wave  (or  waves).  The  understanding  of  backscattering 
properties  may  be  helpful  in  the  evaluation  of  optical  technologies  to  detect 
persistent  bubbles  in  wakes  which  may  be  feasible  in  ocean  regions  (e.g.  the 
Arctic)  where  background  light  scattering  is  less  significant.  (Bubbles  in 
wakes  have  been  asserted  to  have  anomalously  long  lifetimes.)  The 
understanding  of  glory-enhanced  near-forward  scattering  is  related  to  the 
degradation  of  the  polarization  of  circularly  polarized  light  incident  on 
bubbly  liquids. 

The  motivation  for  other  optical  scattering  research  is  noted  in  Sec. 

IVF-H. 


B.  Critical  angle  scattering  of  light  from  bubbles 

Our  principal  recent  work  on  critical-angle  scattering  is  the  analysis 
A4  A10 

and  publication  ’  of  irradiance  measurements  for  the  angular  region  68  to 
85  degrees  and  bubbles  ranging  in  radius  from  26  to  986  pm.  This  angular 
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Fig.  9.  Normalized  scattered  intensity  from  a  bubble  with  ka  =  1633  and  the 
electric  field  parallel  to  the  scattering  plane,  (a)  Three  models:  the  solid 
curve  is  the  Mie  theory,  the  dashed  curve  is  a  physical-optics  approximation, 
and  the  dotted  curve  is  from  simple  geometric  optics,  (b)  The  solid  curve  is 
experimental  data  and  the  dashed  curve  is  the  physical-optics  model. 
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Fig.  10.  The  scattered  intensity  for  ka  612  when  the  electric  field  is 
perpendicular  to  the  scattering  plane.  The  solid  curve  is  data,  the  dotted 
curve  is  the  Mie  result,  and  the  dashed  curve  is  the  physical-optics  model. 
The  data  indicate  the  bubble  radius  was  0.0461  +  0.0001  mm. 


region  includes  the  transition  to  near-total  reflection  from  the  bubble’s 

surface  as  well  as  the  coarse  oscillations  of  the  intensity  which  result  from 

23 

interference  and  diffraction.  The  incident  light  was  a  linearly  polarized 

expanded  beam  of  a  He-Ne  laser.  Figure  9(a)  compares  various  theories  for  the 

scattering  in  this  region  for  a  bubble  having  a  radius  a  =  0.1235  mm  such  that 

ka  =  2iTaAo  =  1633.  The  dashed  curve  is  the  Marston-Kingsbury 
23 

approximation.  The  dotted  curve,  which  is  based  on  geometrical  optics, 

has  an  unphysically-divergent  derivative  at  the  critical  scattering  angle  of 

82.7  deg.  Figure  9(b)  compares  our  measurements  for  a  bubble  having 

a  =  0.12  ram  with  the  results  of  the  Marston-Kingsbury  approximation.  The 

agreement  is  representative  of  measurement  versus  model  comparisons  for  six 

bubbles  in  the  size  range  0.046  mm  to  0.733  mm  published  in  Ref.  A4.  It  is 

also  representative  of  comparisons  based  on  measurements  for  12  other  bubbles 

which  are  displayed  in  Langley's  Ph.D.  dissertation.01  These  comparisons, 

while  confirming  the  basic  features  of  the  Marston-Kingsbury  model,  also 

A4 

reveal  its  nonuniform  nature  for  scattering  angles  >  82  deg. 

For  the  bubble  having  a  =  0.046  mm  measurements  were  compared  with  the 
detailed  results  of  Mie  theory.  This  comparison.  Figure  10,  shows  that  Mie 
theory  correctly  predicts  the  locations  of  both  the  fine  structure  and  the 
contrast  modulations  of  the  fine  structure.  The  success  of  Mie  theory  to 
predict  these  details  of  the  measurements  is  evidence  of  the  sphericity  of  the 
freely  rising  bubble.  The  comparison  also  suggests  that  the  bubble's  radius 
a  was  known  to  within  +0.1  pm.  See  the  Applied  Optics  publication.^ 

One  of  the  difficult  aspects  of  Langley' s  experiment  was  the 
development  of  a  precise  angle  calibration  procedure.  Figure  11  shows  the 
apparatus  during  the  angle  calibration  phase.  Another  important  technique 
used  was  the  trapping  and  direct  radius  measurement  of  single  bubbles  after 
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Fig.  11(b) «  Photograph  of  the  experimental  apparatus.  The  numbered  devices 
are:.  (1)  front  end  of  laser,  (2)  beam  expander,  (3)  polarization  rotator,  (4) 
mirror,  (5)  beam  splitter,  (6)  retroref lecting  prism,  (7)  lens,  (8)  scattering 
chamber,  (9)  camera,  (10)  mechanism  of  angle-calibration  goniometer  and 
rotatable  reflector. 
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PENDANT  BUBBLE  RADIUS  ap  (mm) 


12.  Measurements  and  model  for  the  angular  frequency  of  fine— structure 
lines.  Data  are  displayed  for  120  rising  and  23  pendant  bubbles,  whose  radii 
a  and  a  ,  respectively,  were  measured  by  microscope.  Mie  results  were 
obtained  from  high- re solution  computations. 
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they  had  risen  through  the  scattering  chamber.  These,  and  other  experimental 

A.A 

details,  are  discussed  in  the  Applied  Optics  paper. 

In  addition  to  studying  the  coarse  structure,  Langley  also  measured  the 

angular  spacing  o  the  rapid  oscillations  evident  in  Fig.  9(b)  which  we  call 

the  "fine  structure."  Re  refined  a  model  for  the  fine  structure  quasi- period 

20 

proposed  initially  by  Mars ton.  The  results  of  Langley f s  analysis  are 

compared  with  the  data  in  Fig.  12.  The  comparison  was  carried  out  for  a 

region  centered  on  a  scattering  angle  of  82.7  deg.  The  data  confirms  the 

model  for  both  states  of  incident  polarization  (E— field  parallel-to  or 

perpendicular  to  the  scattering  plane) •  The  smallest  bubble  in  the  rising- 

bubble  data  set  had  a  radius  of  0.026  mm  while  the  largest  bubble  had  a  radius 

of  0.99  mm.  The  comparison  also  includes  data  for  pendant  bubbles  which  were 

attached  to  a  needle.  See  Ref.  A4  for  details. 

We  also  published  a  paper  which  shows  and  discussed  colored  bands  that 

appear  near  the  critical  scattering  angle  when  a  bubble  is  illuminated  with 
A2 

white  light.  For  reasons  of  convenience  the  actual  experiments  were  carried 

out  on  bubbles  in  glass;  however,  the  phenomena  observed  should  be  quite 

similar  to  those  present  for  bubbles  in  water.  Applications  to  the  detection 

and  sizing  of  bubbles  in  water  illuminated  by  white  light  (instead  of  the 

laser  light  used  by  Langley)  are  noted  as  are  applications  to  the  visual 

appearance  of  clouds  of  bubbles.  A  white- light  bubble-size  spectrometer 

12 

proposed  by  MacIntyre  would  make  use  of  a  related  manifestation  of  the 
colored  bands  to  measure  the  population  spectra  of  microbubbles  in  seawater. 

C.  Strong  cross- polarization  in  the  near-forward  scattering  from 
bubbles:  Observations  of  the  forward  optical  glory. 

A  significant  experiment  was  carried  out  which  is  discussed  in  detail 
in  Langley1 s  Ph.D.  dissertation.01  This  was  the  first  observation  of  the 
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Fig.  13.  Apparatus  for  observing  forward  glory  scattering  effects. 
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Fig.  14.  Forward  glory  rays  of  a  bubble  and  associated  focal  parameters. 
Dashed  lines  show  the  origin  of  the  toroidal  wavefront. 
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cross-polarized  forward  optical  glories  of  bubbles  in  water  and  in  silicone 

oil.  (These  experiments  were  done  in  January  -  May  of  1984  and  were  completed 

by  Langley  while  Mars ton  was  away  on  a  sabbatical.)  A  simplified  diagram  of 

the  apparatus  used  is  shown  in  Fig.  13.  The  polarizers  were  crossed  so  that 

in  the  absence  of  a  bubble  a  negligible  amount  of  light  was  transmitted  by 

polarizer  2.  Our  theoretical  analysis^0 >A^’01  shows  that  the  principal  source 

of  near-forward  cross  polarized  scattering  should  be  that  due  to  axially- 

focused  (i.e.  glory)  waves.  The  ordinary  forward-diffraction  peak  which 

contributes  to  the  total  scattering  of  bubbles  (and  also  to  opaque  spheres  or 

discs)  is  intrinsically  co-polarized  with  the  incident  beam  so  that  its 

contribution  to  the  cross— polarized  scattering  may  depend  on  measurement 
A3 

configuration.  Figure  14  illustrates  one  class  of  rays  producing  a  forward- 
facing  axially-focused  toroidal  wave.  (The  reader  may  wish  to  compare  this 
with  Figs.  2  and  8  which  illustrate  rays  leading  to  some  backward  facing 
toroidal  waves.)  In  addition  to  the  2— chord  wavefront  shown  the  lower  3— chord 
ray  also  produces  a  toroidal  wavefront.  Calculations  show  that  the  dominant 
cross-polarized  near-forward  toroidal  wave  is  the  2-chord  wave  for  spherical 
bubbles  either  in  water  (refractive  index  =  1.33)  or  in  silicone  oil 
(refractive  index  =  1.403). 

Most  of  the  data  was  taken  for  bubbles  in  silicone  oil.  This  is 
because  the  high  viscosity  of  the  oil  reduces  to  rise  velocity  of  bubbles  so 
that  it  was  easy  to  work  with  a  single  bubble.  Calculations01  show  that  the 
cross— polarized  irradiance  from  a  bubble  in  water  is  similar  in  magnitude  to 
that  of  a  bubble  in  silicone  oil  and  that  the  general  angular  structure  of  the 
scattering  should  be  similar  apart  from  small  changes  in  scale.  The  far- field 
cross-polarized  irradiance  was  photographed  by  adjusting  the  camera  lens  in 
Fig.  13  such  that  the  film  plane  was  focused  at  infinity.  Figure  15  is  a 
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Wf;17.  Cross-polarized  near- forward  scattering  from  a  cloud  of  small 
bubbles  rising  water. 
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representative  photograph  of  the  scattering  for  a  bubble  in  silicone  oil  for  a 

bubble  whose  radius  a  was  measured  to  be  0.203  mm.  The  quasi- periodic 

nature  of  the  Maltese  cross— like  pattern  is  characteristic  of  cross-polarized 

22  A9 

glories  dominated  by  a  single  toroidal  wave.  *  The  center  of  symmetry 

should  correspond  to  the  direction  of  the  incident  beam.  Except  close  to  the 
forward  direction,  the  angular  spacing  of  the  bright  (or  dark)  rings  should 


Ojp)  radians  ,  (14) 

where  the  ratio  b  /a  =  0.65  is  independent  of  the  bubble  radius  a  and  b 

^  2 

is  the  radius  of  the  focal  circle  as  shown  in  Fig.  14.  Figure  16  compares  the 
measurements  with  the  predictions  of  Eq.  (14).  This  comparison  demonstrates 
that  the  cross-polarized  scattering,  as  shown  in  Fig.  15,  is  dominated  by  the 
2-chord  toroidal  wave  and  not  by  forward  diffraction.  Forward  diffraction 
should  produce  angular  structures  having  a  quasi-period  of 

forward  diffraction  ka  ra<^ans  *  (15) 

which  clearly  differs  from  that  of  Eq.  (14)  and  of  the. data.  (Structures 
having  a  quasi-period  given  by  Eq.  (15)  should,  however,  dominate  the  co¬ 
polarized  scattering.) 

For  bubbles  in  water,  it  was  not  feasible  to  photograph  the  scattering 
from  a  single  bubble  and  then  to  trap  it.  Instead  the  cross- polarized 
scattering  of  a  cloud  of  small  freely-rising  bubbles  in  water  was 
photographed.  Figure  17  is  a  representative  photograph  of  the  scattering. 

For  bubbles  in  water  the  quasi-period  should  be  given  by  Eq.  (14)  with 
a  =  0*689*  In  Fig.  17,  quasi— periodic  oscillations  are  evident  with 
A<j>  -  0.005  radian  -  0.29  deg.  From  Eq.  (14)  we  estimate  the  average  bubble 
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(a) 


Figure  18 
Photographs  of  the 

(a)  co-polarized  and 

(b)  cross-polarized 
forward  glory  circles 
for  a  bubble  with  a 
radius  of  1.88  am  in 

a  liquid  of  refractive 
index  1.403.. 


(b) 


radius  of  the  cloud  to  be  ^  69  pm.  These  observations  indicate  that  bubbles 
of  this  radius  are  sufficiently  round  to  manifest  a  forward-optical  glory. 
Bubbles  in  water  become  oblate  if  the  radius  is  somewhat  larger  than  this 
value  so  it  is  not  known  how  large  a  bubble  may  be  and  still  manifest  the 
glory. 

Cross-polarized  optical  glories  should  also  be  present  when  the 
incident  light  is  circularly  polarized.  They  should  appear,  as  in  Fig.  15  and 
17  except  that  the  irradiance  will  not  depend  on  the  azimuthal  angle  so  that 
the  Maltese  cross- like  pattern  will  be  lost.  The  pattern  should  be  quasi¬ 
period  in  scattering  angle.  The  relation  between  scattering  patterns  with 
linearly  and  circularly  polarized  incident  light  are  derived  by  us  in  Ref.  A9 
for  near  backward  scattering.  We  have  also  derived  similar  ones  for  near- 
forward  scattering. 

Langley  also  photographed01  the  focal  circles  associated  with  the 
forward  facing  glory  waves.  This  was  done  by  focusing  the  camera  in  Fig.  13 
directly  on  the  bubble  so  that  the  focal  circles  appear  as  thin  bright  rings. 
The  bubbles  were  in  silicone  oil.  Figure  18  shows  enlargements  of 
representative  photographs.  In  each  case,  more  than  one  focal  circle  is 
visible  corresponding  to  the  sources  of  glory  waves  having  different  numbers 
of  chords  as  in  Fig.  14.  The  ratios  of  the  ring  radii  to  that  of  the  bubble 
radii  were  measured  and  were  found  to  agree  with  predictions.01 

D.  Strong  cross-polarization  in  the  near- backward  scattering  from 
bubbles  in  water:  The  backward  optical  glory. 

The  principal  advancements  in  this  aspect  of  the  research  during  this 
contract  period  were  improvements  to  the  theory.  Observations  were  previously 
made  by  us  (with  the  support  of  this  contract)  which  confirmed  the  existence 
of  backward  glory  scattering  from  bubbles  in  silicone  oil22  and  in  water. A9 


y  (degrees) 

Fig.  19.  Theoretical  crosa-polari zed  (1^^)  and  co-polarized  (1^'^)  optical 
irradiances  plotted  as  a  function  of  the  backscattering  angle  y  for  a  bubble 
in  water  of  radius  -  66  pm  (for  green  light).  The  solid  and  dashed  curves  are 
exact  for  a  spherical  bubble  and  were  calculated  from  Mie  theory.  The  dotted 
and  dot-dashed  curves  are  from  a  physical-optics  approximation  of  the  axially- 
focused  (glory)  waves  and  the  specular  axial  reflections.  Some  of  the  rays 
included  in  the  calculation  are  shown  in  the  inset. 
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For  bubbles  in  silicone  oil  the  cross-polarized  patterns  were  like  that  of 
Fig.  15. 

The  first  theoretical  advancement  was  to  formulate  a  novel  uniform 

analysis  of  cross-polarized  and  co-polarized  irradiances  according  to  Mie 
A3 

theory.  This  short  paper  concerns  analytical  details  and  definitions,  which 

though  important,  need  not  be  repeated  here. 

The  other  theoretical  advancement  is  more  physical  in  nature  and  merits 

description.  For  bubbles  in  silicone  oil  the  cross-polarized  backward  glory 

22 

is  dominated  by  the  3— chord  glory  wave.  For  bubbles  in  water,  it  was 

evident,  however,  that  it  would  be  necessary  to  sum  the  electric  fields  due  to 

several  different  glory  waves  in  order  to  describe  the  near-backward 

scattering.  Such  a  summation  was  now  possible  because  of  the  analysis  of  the 

angle-dependent  phase  shift  of  glory  waves  first  carried  out  by  Marston  and 

described  in  the  paper  on  the  acoustic  glory  of  fluid  spheres. ^  The  resulting 

physical-optics  model  predicts  irradiances  which  agree  with  those  given  by  Mie 

A9 

theory  for  spherical  bubbles  in  water.  This  confirms  several  details  of  the 
physical— optics  model  of  the  optical  glory  of  bubbles.  Figure  19  compares  Mie 
theory  (solid  and  dashed  curves)  with  model  results  (dotted  and  dot-dashed 
curves) .  The  scattering  is  expressed  in  normalized  units  in  which  I  is  the 
irradiance  which  would  be  reflected  by  a  perfectly  reflecting  sphere  of  the 
same  size.  The  solid  and  dashed  curves  apply  to  cross-polarized  and  co¬ 
polarized  irradiances  and  they  are  to  be  used  with  the  right  and  left  vertical 
scales,  respectively.  The  incident  light  is  linearly  polarized  and  the 
calculation  is  for  an  aximuthal  angle  of  45  degrees  relative  to  its  direction 
of  polarization.  For  ka  =  1000,  a  bubble  in  water  illuminated  by  green  light 
has  a  radius  of  about  66  ym.  The  inset  on  the  upper  right  of  Fig.  19  shows 
some  of  the  important  rays . 
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E.  Use  of  light  scattering  to  measure  the  amplitude  of  radial 

oscillations  of  bubbles  in  water:  An  investigation  of  the  theory. 

In  spite  of  considerable  theoretical  interest  and  practical 

applications,  direct  measurements  of  the  response  of  small  buobles  in  water  to 

sound  of  moderate  or  low  amplitudes  have  not  been  made.  Such  measurements 

would  be  important  since  it  has  been  recently  realized  that  the  response  may 

12 

be  affected  by  the  presence  of  surface  monolayers  and  by  evaporation  and 
condensation,  both  of  which  are  difficult  to  model.  A  direct  measurement 
would  require  an  optical  technique  which  was  sensitive  to  radial  pulsations 
with  amplitudes  Aa  of  <  2  |im.  A  technique  with  the  required  sensitivity  was 
suggested  by  Marston  in  his  review  presentation  of  light  scattering  from 
bubbles  given  at  the  Workshop  on  Wave  Propagation  in  Bubbly  Liquids  (Rosslyn, 
VA,  May  1983)  which  was  sponsored  by  the  Office  of  Naval  Research.  The 
proposed  technique,  which  would  use  features  of  the  coarse  structure  in  near 
critical-angle  scattering  (see  Sec.  IVB) ,  is  outlined  below.38 

Stimulated  in  part  by  our  descriptions  ’  *  of  scattering  from 

bubbles,  G.  M.  Hansen  and  Professor  L.  A.  Crum  of  the  University  of 
Mississippi  recently  used  the  coarse  structure  to  calibrate  a  photometer  and 
successfully  measure  the  average  radius  of  bubbles  in  an  acoustic  levitator.2^ 
Thev  are  also  interested  in  direct  light  scattering  measurements  of  a  bubble's 
oscillation  amplitude .  Consequently  the  theoretical  comparison  of  techniques 
outlined  below  is  merited. 

Perhaps  the  simplest  technique  for  detecting  small  oscillations  would 
be  to  scatter  light  at  an  angle  close  to,  or  equal  to,  the  critical  scattering 
angle  of  82.7  deg.  The  Marston-Kingsbury  approximation23  predicts  that  the 
flux  into  a  photometer  which  collects  the  light  in  a  small  angular  region 
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containing  this  angle  will  be  proportional  to  the  critical-angle  irradiance 
I 


c 


IC  ~  4  XR  s  IR  Iinc(2R^  * 


(16a, b) 


where  I  is  the  irradiance  scattered  by  a  perfectly  reflecting  sphere  of  the 
same  radius  a,  1^  is  the  incident  irradiance,  and  R  is  the  distance  to 
the  photometer  from  the  bubble.  (Corrections  due  to  windows  are  assumed  to  be 
time- independent  and  will  be  neglected.)  In  response  to  a  small  change  in  a 
by  a,  the  resulting  shift  in  irradiance  is  due  to  the  change  in  geometric 
cross-section.  From  Eq.  (16)  it  is  evidently  given  by 


Alc  -  i  IR  (Aa) / a  , 


(17) 


since  dl  /da  =  2  1,,/a. 

K.  R 

A  more  sensitive  (but  more  difficult)  method  of  measuring  a  would  be 
to  place  the  photometer  at  an  angle  somewhat  less  than  82.7  deg.  The  observed 
region  should  be  close  to  the  center  of  the  Nth  course  oscillation  below  82.7 
deg  where  the  slope  of  the  irradiance— versus— scattering— angle  curve  is 
predicted  to  be  large.  For  example  one  might  use  70  deg  for  the  bubble  in 
Fig.  9.  The  irradiance  change  will  be  enhanced  since  small  changes  in  the 
radius  will  shift  the  course  structure.  When  a  increases  by  Aa,  the  coarse 
feature  observed  has  an  angle  shift  A(p  of  approximately 


(18) 


1 

20  23  —~o~ 

where  we  use  the  result  *  that  the  quasi- period  ft  is  roughly  *  a~  .  Let 


I  denote  the  peak-to-peak  amplitude  of  the  Nth  coarse  cycle*  In  response 
to  the  change  Aa,  the  irradiance  change  AI  is  roughly 
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AI  -  2tt 


I  A<f> 
PP 
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(19) 


For  the  aforementioned  example  from  Fig.  9  the  Mie  theory  curve  gives 

I  -0.7  I,,  and  the  photometer  was  assumed  to  be  set  where  N  =  2.  Then  a 
PP  R 

comparison  of  Eq.  (19)  and  (17)  gives  AI  -  4.4  AI  •  Additional  sensitivity 

c 

could  be  achieved  by  placing  the  photometer  at  even  smaller  scattering  angles 
since  both  N  and  I  would  be  larger*  Furthermore,  it  may  be  argued  that 
Eq.  (19)  underestimates  AI  since  the  increase  in  the  geometric  cross-section 

should  give  an  additional  term  in  Eq.  (19)  of  magnitude  2  1^  (Aa)/a.  In  the 

.  R 

final  analysis,  It  might  be  necessary  to  calibrate  the  photometer  using 
calculated  Mie  Intensities  averaged  over  the  aperture  of  the  photometer.  It 
would  be  necessary  for  experimentalists  to  measure  both  the  constant  and  the 
oscillating  components  of  the  photometer's  output.  N  may  be  determined  by 
counting  the  coarse  variations  in  the  output  as  the  angle  was  varied  from  82.7 
deg.  to  the  final  value. 

F.  Rainbow-enhanced  backscattering  from  a  glass  sphere  in  water: 

Theory  of  a  novel  aspect- insensitive  retro- ref lector  for  use 
in  water. 

In  our  paper  on  the  acoustical  glory  from  fluid  spheres^*  a  novel  method 
for  enhanced  scattering  was  proposed  which  combined  axial  focusing  with  the 
weak  focusing  of  rainbow  or  "stationary”  ray.  It  was  also  noted  there,  that 
this  dual  focusing  should  enhance  the  light  backscattered  from  dielectric 
spheres  whose  relative  optical  refractive  indices  were  about  1.180  and  1.089. 
These  relative  indices  would  be  realizable  by  placing  glass  spheres  in  water. 

A  physical-optics  model  of  this  dual  focusing,  which  also  describes  the 


depolarization  of  the  scattered  electromagnetic  wave,  was  subsequently 

formulated  and  confirmed  by  comparison  with  the  results  of  exact  Mie  theory. 

A9 

This  research  was  summarized  in  the  manuscript  for  the  Ocean  Optics 
VII  conference.  The  polarization  properties  and  (large)  magnitude  of  the 
scattering  illustrated  by  figures  in  that  paper.  Applications  to  the  design 
of  optical  tracers,  which  could  be  made  neutrally  buoyant  in  sea  water,  are 
discussed 

G.  Novel  method  for  the  discovery  and  classification  of  focused 

waves:  Hyperbolic-umbilic  diffraction  catastrophe  in  the  light 
scattered  by  an  acoustically-levitated  drop. 

qo 

There  has  been  recent  progress  in  the  calculation  of  optical  and 

acoustical10  scattering  patterns  from  penetrable  spheroids  at  long 

wavelengths.  At  short  wavelengths,  however,  partial-waves  sums  are  difficult 

to  evaluate  and  the  scattering  amplitudes  may  vary  so  rapidly  with  angle  that 

the  gross  features  may  be  difficult  to  identify.  Furthermore,  even  at  short 

wavelengths,  geometric  optics  Is  not  applicable  near  angular  foci  where  the 

scattering  is  enhanced.  This  is  because  at  foci,  simple  stationary  phase 

33 

approximations  of  diffraction  integrals  diverge.  Consequently,  the 

discovery  and  classification  of  angular  foci  is  important  for  the 

understanding  of  acoustical  and  optical  short-wavelength  scattering.  The 

theoretical  framework  for  the  classification  of  foci  was  established  recently 

33  3 A 

(mostly  by  European'  and  Soviet  researchers)  and  is  known  as  "catastrophe 
optics . ” 

In  a  recent  experiment  Marston  discovered  a  new  focus  or  "diffraction 

catastrophe"  in  the  light  scattered  by  an  acoustically  levitated  oblate  drop 
A12  B9 

of  water  in  air.  *  A  concise  description  of  the  experiment  (a  preprint  of 
reference  A12)  is  included  in  the  present  report  as  Sec.  VIII  the  Appendix. 


In  addition  to  the  results  discussed  there,  the  significance  of  this  type  of 

experiment  is  as  follows:  (i)  Once  the  scattering  patterns  (Appendix,  Fig.  4) 

were  observed,  it  was  evident  that  the  focus  must  be  a  hyperbolic-umbilic 

diffraction  catastrophe.  It  follows  from  previous  research  in  catastrophe 

optics  that  no  other  class  of  focusing  could  explain  these  patterns,  (ii)  It 

is  not  surprising  that  our  experiments  preceded  our  consideration  of  the 

theory  since  ray  tracing  through  (or  ray  reflection  from)  complicated  shapes 

is  difficult.  It  is  often  difficult  to  predict  the  location  and  types  of  foci 

in  advance  of  a  difficult  calculation.  Consequently,  experimental  searches 

for  diffraction  catastrophes  in  scattering  patterns  are  worthwhile,  (iii) 

Optical  analogs  may  be  used  in  place  of  acoustical  scatterers.  For  penetrable 

objects,  it  is  of  course  necessary  that  the  optical  refractive  index 

correspond  to  the  acoustic  refractive  index  of  interest,  (iv)  Once  the 

catastrophe  type  has  been  classified,  scaling  laws  may  be  used  to  predict 

the  scattering  for  other  short  wavelengths.  For  example,  the  intensity  at  the 

hyperbolic-umbilic  focus  (Fig.  4d  of  the  Appendix)  must  be  proportional  to 
2/3  2 

(D/A)  (D/r)  where  D  is  the  drop1  s  diameter,  A  is  the  wavelength,  and 

r  is  the  distance  from  the  drop. 

This  experiment  was  carried  out  with  the  collaboration  of  Dr.  E.  Trinh 
while  Marston  was  on  a  sabbatical  leave  (from  April- June  1984)  at  the  Jet 
Propulsion  Laboratory  (Pasadena,  CA) .  Analysis  of  that  data  was  performed  by 
Marston  after  returning  to  Washington  State  University  under  the  support  of 
the  present  ONR  Contract. 


H.  Radiation  torque  on  a  sphere  caused  by  a  circularly-polarized 
electromagnetic  wave. 

As  previously  noted  in  Sec.  IVB  we  have  calculated  the  scattering  of 
circularly  polarized  light  from  bubbles/9  During  the  course  of  that 
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research,  the  following  question  was  evident:  "Would  it  be  possible  to  use 
the  angular  momentum  transport  of  photons  (in  e.g.  a  circularly-polarized 
laser  beam)  to  spin-up  a  sphere?"  Fundamental  to  this  question  is  the 
calculation  of  the  radiation  torque  on  a  sphere  (of  arbitrary  size)  caused  by 
a  circularly-polarized  electromagnetic  wave.  This  fundamental  calculation  had 
not  been  previously  published  even  though  there  have  been  several 
investigations  of  the  radiation  pressure  on  spheres.  Some  potential 
applications  of  such  a  calculation  include:  (i)  The  evaluation  of  proposed 
methods  for  "burning"  or  "clearing”  holes  through  fog  with  an  intense  laser 
beam.  (The  transport  of  angular  momentum  to  drops  should  enhance  the  clearing 
of  holes.)  (ii)  The  estimate  of  instrumental  bias  in  a  proposed  method  for 
measuring  fluid  vorticity.  (The  method,  based  on  the  scattering  of  circularly 
polarized  light  from  small  spheres,  may  be  useful  in  acoustics  and 
hydrodynamics  research  or  in  remote  sensing.)  (iii)  A  novel  bolometer  for 
circularly  polarized  microwave  radiation. 

A  direct  calculation  of  the  torque  was  carried  out  by  Marston  and 
Professor  J.  H.  Crichton  of  Seattle  Pacific  University.  A  summary  has  been 
published^  and  the  detailed  calculation  has  been  accepted  for  publication  in 
Physical  Review. ^  The  simplest  result  is  for  a  z-directed  plane  wave  with 
pure  circular  polarization  corresponding  to  a  positive  value  of  the  photon 
spin.  Formulation  of  the  Maxwell  stress  dyad  of  the  total  (incident  and 
scattered)  field  gives  the  following  torque  relative  to  the  sphereTs  center 


rz  -  Ira2W“ 


Here  I  and  oo  are  the  incident  wave*  s  irradiance  and  angular  frequency  and 

a  and  0  ,  are  the  sphere's  radius  and  Mie- theoretic  absorption  efficiency. 
xabs 

Consequently  the  effective  cross-section  for  torque  is  the  same  as  that  for 
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energy  absorption  Tra  Qa^s*  The  Pro°f  is  complicated  since  it  involves 

detailed  consideration  of  the  angular  momentum  transport  of  multipole  fields 

and  the  correction  of  some  previous  scientific  literature*  We  have  also 

estimated  the  steady-state  angular  velocity  achieved  when  the  torque  is 

A7 

balanced  by  viscous  drag  and  the  time  required  to  spin-up  to  this  angular 

i  ^  B6 
velocity* 

Since  the  torque  was  found  to  be  linked  to  the  absorption  of  energy,  it 
now  appears  that  the  effect  may  be  too  small  to  enhance  the  clearing  of 
optical  holes  through  fog*  The  absorption  of  radiation  by  the  drop  is  more 
likely  to  be  significant  to  the  clearing  process  than  is  the  torque  which 
accompanies  it;  however,  a  detailed  analysis  has  not  been  performed. 

V.  Characterization  of  Photo-Acoustic  Sources 

A.  Production  of  sound  by  a  pre-existent  bubble  in  water  illuminated 

by  modulated  light:  A  novel  photo- acoustic  source 

One  of  the  goals  of  this  research  is  to  examine  whether  pre-existent 

bubbles  can  enhance  the  conversion  of  modulated  optical  energy  to  acoustical 

energy*  Our  proposed  mechanism  for  this  energy  conversion  are:  (i)  the 

compression  of  the  bubbles  by  optical  radiation  pressure  and  the  resulting 

radiation  of  sound  from  the  monopole  oscillations  of  the  bubble;  (ii) 

radiation  associated  with  monopole  oscillations  driven  by  the  heating  and 

cooling  of  the  gas  within  the  bubble*  (The  heating  being  due  to  the 

absorption  of  light  by  the  gas  and  by  the  bubblers  surface  "skin"  or 

surroundings »)  We  have  estimated  that  the  relevant  radiation  pressure  is 

primarily  due  to  light  reflected  from  the  bubble1 s  surface  for  which  the  angle 

of  incidence  exceeds  the  critical  value  0  »  arcsin(n  -i)  where  n  is  the 

c  o  o 

liquid’s  refractive  index.  When  averaged  over  the  bubble's  surface,  the 
radial  projection  of  the  radiation  pressure  becomes 
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I  3q 

p  -  —  n  cos  0  , 

r  3c  o  c  ’ 


(20) 


where  I  is  the  optical  irradiance  (watt/ta  )  and  c  is  the  speed  of  light 
in  a  vacuum.  To  estimate  the  maximum  radiated  pressure,  assume  that  the 
irradiance  is  modulated  in  such  a  way  that 


pr  =  pro  ^  +  cos  » 


(21) 


and  choose  u)  to  correspond  to  the  monopole  resonance  frequency  of  the 
bubble.  The  amplitude  of  the  sound  wave  a  distance  r  from  a  bubble  of 
radius  a  is 


ps  *  pr0  0  r  ’ 


(22) 


where  Q  is  the  quality  factor  of  the  resonance.  For  a  bubble  in  water  having 

a  =  20  urn,  Q  -  10  and  we  anticipate  |p  |  =  10- 1  Pa  for  our  apparatus.  This 

s 

amplitude,  though  small,  is  not  smaller  than  signals  detected  in 
35 

conventional  photo-acoustic  experiments.  It  should  be  emphasized  that 
radiation  due  to  modulated  heating  of  the  bubble  may  greatly  exceed  that 
predicted  by  (20)  and  (22)  due  to  radiation  pressure.  The  derivations  of  Eqs. 
(20)  and  (22)  were  carried  out  by  Marston  and  Unger  during  a  previous  contract 
period.  The  fundamental  photoacoustic  mechanism  proposed  is  different  than 
previously  known  mechanisms  which  are  reviewed  in  Ref  35 . 

To  facilitate  the  detection  of  small  signals  radiated  by  single 
bubbles,  it  is  necessary  to  "levitate"  single  bubbles  in  a  quiet  environment. 
In  1983  we  demonstrated  a  technique  for  levitating  bubbles  by  counteracting 
their  buoyancy  with  the  radiation  force  of  a  downward  propagating  laser  beam. 
It  can  be  shown  that  a  stable  levitation  position  can  exist  when  laser  is 
operated  in  the  TEM*01  mode  which  is  usually  described  as  the  "doughnut"  mode 
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Fig.  20.  Experimental  evidence  of  the  radiation  of  sound  from  a  gas  bubble 
(in  water)  subjected  to  a  burst  of  light  pulses  of  peak  power  =  1  watt:  (a) 
shows  the  electrical  input  to  the  optical  modulator  consisting  of  3  equally 
spaced  pulses  each  of  duration  =2.6  ysec;  (b)  shows  the  hydrophone  signal 
output  when  a  bubble  is  present;  and  (c)  shows  the  hydrophone  signal  output 
when  no  bubble  is  present.  The  trace  in  each  case  is  synchronized  to  2.6  ysec 
prior  to  the  start  of  the  electrical  input  as  shown  in  (a).  Figure  21  shows 
the  apparatus . 


Laser 


Fig.  21.  Diagram  of  the  appara tue  used  both  to  optically  levitate  bubbles  and 
to  subject  them  to  bursts  of  light  pul sea#  The  receiver  hydrophone  tone lata 
of  a  ring  of  thin  piezoelectric  material  concentric  with  the  light  beam. 
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because  of  an  intensity  minimum  in  the  center  of  the  beam*  With  the  output  of 
a  3  watt  cw  Ar-Ion  laser  (purchased  with  the  combined  support  of  contract  and 
Washington  State  University  funds) ,  Unger  has  levitated  several  bubbles  in  the 
size  range  with  radii  of  10  to  25  ym  for  periods  up  to,  and  occasionally 
exceeding,  15  minutes.  The  bubbles  contain  hydrogen  (H^)  gas  and  are 
generated  in  the  water  by  electrolysis  at  the  surface  of  a  thin  tungsten  wire. 

A  simplified  diagram  of  the  apparatus  for  detecting  the  photo-acoustic 
response  is  shown  in  Fig,  21.  The  beam-splitter  cube  BS  allows  the  beams 
from  two  lasers  to  be  combined  without  significant  power  losses  (because  of 
the  polarization  properties  of  the  cube  and  of  the  lasers)  .  The  combined 
beams  are  coaxial.  Laser  No.  1  is  the  aforementioned  Ar-Ion  laser  and  its 
beam  is  used  to  trap  and  levitate  a  single  gas  bubble  in  the  water  tank. 

Laser  No.  2  is  a  second  Ar-Ion  laser  made  available  with  contract  funds.  Its 
cw  beam  (about  1  watt)  is  modulated  using  an  electro-optic  modulator  system 
(Conoptics  Model  10,  Danbury,  CT;  purchased  with  contract  funds).  The 
modulated  output  is  usually  taken  to  be  a  burst  of  light  pulses  similar  In 
form  to  the  eight-pulse  burst  shown  in  Fig.  20(a).  (Actually,  Fig.  20(a) 
shows  the  electrical  input  to  the  modulator.  The  modulator  is  adjusted  so 
that  the  intensity  history  of  the  burst,  monitored  with  a  photo  tube  is 
similar  to  the  input  shown.)  The  frequency  of  pulses  in  the  burst  in  this 
case  was  190  kHz.  The  bandwidth  of  the  modulator  is  DC  to  10  MHz.  The  reason 
for  using  a  short  burst  is  so  that  acoustic  reflections  from  the  surfaces  of 
the  tank  would  not  interfere  with  the  direct  sound  from  the  bubble.  The  burst 
was  repeated  3  times  per  second. 

Figure  20(b)  shows  preliminary  results  which  appear  to  demonstrate  a 
bubble  acts  as  a  photoacoustic  source.  It  shows  the  output  of  the 
preamplifier  averaged  over  256  repetitions  of  the  input  burst  shown  In  Fig. 


20(a).  That  the  signal  in  Fig.  20(b)  was  associated  with  the  presence  of  a 

levitated  bubble  was  confirmed:  Fig.  20(c)  shows  an  average  obtained  wit  the 

same  procedure,  but  immediately  after  the  bubble  had  drifted  out  of  the 

levitation  trap.  Figure  20(c)  shows  no  signal  of  comparable  magnitude  to  the 

one  in  Fig.  20(b).  This  is  to  be  expected  since  their  was  negligible 

absorption  of  the  green  light  used  in  this  experiment.  The  amplitude  of  the 

sound  generated  by  weak  absorption  and  the  usual  thermal-photo-acoustic 
35-37 

mechanism  is  estimated  to  be  small  in  comparison  to  the  bubble's  signal 

as  estimated  from  Eq.  (21).  Comparison  of  Figs.  20(b)  and  20(a)  shows  that 
signal  delay  is  close  to  the  expected  propagation  delay  of  about  13  psec  for 
travel  from  the  bubble  to  the  hydrophone  (a  distance  =  20  mm  in  water). 

For  technical  reasons,  the  radius  of  the  bubble  levitated  during  the 
acquisition  of  Fig.  20(b)  could  not  be  measured;  however,  from  properties  of 
the  source  for  this  bubble,  we  estimate  it  to  be  about  20  Urn.  The  monopole 
resonance  frequency  of  such  a  bubble  is  about  160  kHz.  The  reason  why  the 
radius  is  not  known  is  that  the  transducer  blocks  the  direct  (side)  view  of 
the  bubble.  Various  methods  for  measuring  the  bubble's  radius  (with  the 
transducer  in  place)  are  being  tried,  such  as  use  of  the  coarse  structure  in 
the  scattering  discussed  in  Sec.  IVB  of  this  report.  The  data  in  Fig.  20(b) 
are  representative  of  the  two  cases  in  which  all  of  the  subsystems 
(levitation,  electro-optic,  electro-optic  modulator,  and  signal  processing) 
have  worked  simultaneously.  Attempts  to  obtain  additional  repetitions  of 
this  measurement  are  being  carried  out  as  of  the  time  of  this  writing. 

B.  Transient  photo- acoustic  pulses  generated  by  absorption  and 
heating  in  water 

The  common  mechanism  for  photo— acoustic  pulse  production  involves 
transient  heating  and  relaxation  of  the  liquid  following  the  bulk  absorption 
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of  a  laser  pulse.  To  be  able  to  detect  a  signal  by  this  mechanism  with  our 
apparatus  (Fig.  21  with  Laser  1  turned  off),  it  was  necessary  to  add  a  blue 
dye  to  the  water.  This  is  because  pure  water  is  nearly  transparent  to  the 
green  laser  light.  Our  experiments  typically  used  single  pulses  (of 
duration  ~  10  usee)  or  bursts  of  pulses  such  as  that  shown  in  Fig.  20(a). 

They  were  principally  carried  out  to  confirm  the  magnitude  of  the  sensitivity 
of  our  hydrophones  and  as  a  precursor  to  detailed  hydrophone  calibration 
experiments ♦ 

To  aid  in  the  characterization  of  these  thermal- photoacoustic  line 

sources,  Marston  examined  the  theory  for  the  pressure  profile  radiated  by  a 

A5 

narrow  light  beam  in  a  weak  optical  absorber.  It  turns  out  (see  Ref.  A5), 
for  conditions  close  to  those  of  our  apparatus  where  the  beam  radius  R  and 
pulse  duration  At  are  such  that 

(4DAt)1/2  <  R  «  cAt  ,  (23) 

the  pressure  at  distances  r  »  cAt  may  be  approximated  by 

P<r.t)  -  (g~) 1/2  (^)1/2  s(t  -  r/c)  ,  (24) 

where  the  half- order  derivative  operator  is  defined  in  Sec.  IIIF,  Eq.  (13a), 
and  the  instantaneous  mass  outflow  per  length-of-source  s(t)  is 

s(t)  =  a3P(t)  Cp"1  ,  (25) 

where  P(t)  is  the  instantaneous  beam  power,  a  is  the  optical  absorption 
(for  which  the  P  =  exp(-az)  where  z  is  the  propagation  distance),  3  is 
the  thermal  coefficient  of  volume  expansion,  and  C  is  the  specific  heat  at 


constant  pressure.  In  Eqs.  (23)  and  (24),  c  is  the  velocity  of  sound  in 
-3  2  -1 

water  and  D  -  10  cm  sec  is  the  thermal  diffusivity. 

1  /2 

The  evaluation  of  (d/dt)  si/2’  where  s\/2  is  defined  as  in  E<?’ 
(13b),  was  given  by  Harston  in  Ref.  A5.  It  would  be  applicable  to  laser 
pulses  having  P(t)  «  s^^;  however,  this  is  not  a  common  time  profile  for  a 
modulated  laser  pulse.  On— the-otherhand ,  intensity  histories  from  electro- 
optic  modulators  may  be  obtained  which  are  well  approximated  by  a 
superposition  of  haversine-wave  bursts  of  the  form 

hOT(t)  =  H(ut)(l  -  coswt)  ,  (26) 

where  H  is  the  Heaviside  unit-step  function  defined  in  Sec.  HIE  and  to  is 
radian  frequency.  Equation  (26)  represents  a  burst  of  half- infinite  duration. 
An  N  cycle  haversine  burst  is  given  by 

\  3  hjt)  -  hOT  (t  -  co"1  2TTN)  ,  (27) 

for  integer  N.  To  predict  the  profile  of  such  a  burst,  the  student  working  on 
this  project  (Unger)  has  calculated  that 

(d/dt)1^  h^  =  <o1//2  H(<ot)[sin(o>t-ir/4)  +  21/2f(u)]  ,  (28) 

1  /2 

where  u  =  (2  (ot/fr)  and  f  is  an  auxiliary  Fresnel  integral  function 

38 

defined  by  Eq.  (7.3.5)  of  Abramowitz  and  Stegun.  The  pressure  profiles 

predicted  for  haversine  bursts  have  smoother  leading  and  trailing  edges  than 

1  /2 

the  profiles  of  (d/dt)  s^  published  in  Ref.  A5. 

A  detailed  experimental  confirmation  of  the  aforementioned  profile 
calculation  has  not  yet  been  possible  because  of  drift  of  the  electro-optic 
modulator’s  fidelity.  (This  problem  is  judged  not  to  be  insurmountable; 
however,  stability  is  important  since  the  hydrophone's  output  must  always  be 
averaged  over  repeated  bursts.)  Note  that  other  recently  published 
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measurements  *  of  photo-acoustic  pressure  profiles  were  for  highr-power 
laser  pulses  of  short  duration  (typically  less  than  0.1  nsec)  .  In  those 
measurements,  the  condition  R  «  c  At  of  Eq.  (23)  did  not  hold  so  that  the 
profiles  may  not  have  been  describable  by  Eqs.  (24)  and  (25). 

VI.  Rapid  Cavitation  Induced  by  the  Reflection  of  Shock  Waves 

Progress  on  this  experiment  subsequent  to  the  publication  of  Ref.  39 
was  described  in  an  Ad-Hoc  Report  issued  January  6,  1984  under  ONR  Contract: 
N00014-80— C— 0838 .  The  level  of  support  for  the  present  contract  (and  our 
efforts  to  make  significant  progress  on  the  research  discussed  in  Sec.  III-V 
of  the  present  report)  were  such  that  there  was  no  significant  progress  on 
this  problem  since  the  Ad-Hoc  Report  was  issued.  The  principal 
accomplishments  noted  in  this  report  are  summarized  below.  Technical  aspects 
are  discussed  in  more  detail  in  the  report  (which  is  available  from  P. 

Marston) . 

The  research  is  summarized  as  follows : 

(i)  A  novel  technique  was  demonstrated  for  the  measurement  of  the 
response  of  water  to  short  duration  (typically  1.7  ys)  pulses  of  tension.  The 
tension  pulses  are  created  by  reflecting  a  shock  (i.e.  compression)  pulse  from 
a  quasifree  water-Mylar-air  interface.  The  shock  pulse  is  created  when  an 
impactor  strikes  a  buffer  plate,  the  water  is  located  behind  the  buffer.  See 
Fig.  22  and  Ref.  39.  The  technical  difficulties  described  in  Ref.  39  (due  to 
the  separation  of  the  rear  Mylar  surface  from  the  water)  appear  to  have  been 
overcome  by  reducing  the  thickness  of  the  rear  Mylar. 

(ii)  The  technique  was  used  to  measure  dynamic  tensile  strengths  of 
moderately  pure  (filtered  and  degassed)  water  and  ethylene  glycol.  The 
majority  of  the  water  shots  indicate  that  there  was  a  significant  amount  of 
cavitation  within  1  ys  after  the  creation  of  tension  (in  the  absence  of 
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.  ig •  22*  Diagrams  of  the  impactor,  target,  and  displacement  optical— 
interferometer  used  in  the  shock  wave  experiment.  Rapid  cavitation  was 
induced  in  the  water  following  the  reflection  of  a  shock  pulse  from  the  water- 
Mylar-air  interface  of  the  target. 
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cavitation,  the  negative  pressures  would  have  approached  -13  MPa).  The 
inferred  dynamic  tensile  strength  of  the  water  was  typically  between  4  and  6 
MPa  (i.e.  40  to  60  atmospheres).  Shot  4-01-82  gave  a  dynamic  tensile  strength 
of  11.5  MPa. 

The  data  for  ethylene  glycol  are  more  varied.  The  inferred  dynamic 
tensile  strengths  are  in  excess  of  4  MPa. 

(iii)  The  impact  technique  developed  gives  meaningful  results  as  there 
appears  to  be  no  other  way  to  interpret  the  majority  of  the  shot  records  than 
that  given  as  item  (ii)  (above).  Furthermore,  there  is  clear  evidence  of 
reverberation  of  the  pressure  wave  radiated  by  the  bubbles  within  the  spalled 
layer.  The  impact  technique,  however,  is  not  easy  to  use.  This  is  because  of 
the  low  value  of  the  impact  velocity  (typically  10  m/s)  required  for  the 
production  of  compression  and  tension  of  relatively  small  amplitude  (typically 
13  MPa) .  By  way  of  comparison,  the  typical  impact  experiment  carried  out  at 
the  Washington  State  University  Shock  Dynamics  Laboratory  has  an  impact 
velocity  ~  1000  m/s  and  a  pressure  ~  2  GPa.  The  difficulties  with  low- 
velocity  impact  experiments  are  principally  two- fold: 

(a)  A  small  tilt  of  the  impactor  produces  a  large  tile  of  the  shock  and 
rarefaction  waves.  Such  a  tilt  obscures  the  quantitative 
interpretation  of  the  experiments.  (An  approximation  for  the 
increase  in  tilt  is  noted  in  footnote  13  of  Ref.  39).  To  obtain 
meaningful  experimental  results,  it  is  necessary  that  the  tilt  of 
the  impactor  be  <  1  mrad  —  0.06".  This  requirement  necessitated 
the  development  of  a  novel  optical  alignment  procedure. 

(b)  Reliable  triggering  of  the  oscilloscopes  for  the  recording  of  data 
during  the  desired  time  window  (which  is  a  few  s  in  duration)  is 
more  difficult  in  low-velocity  shots.  This  difficulty  was  overcome 
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by  use  of  two  independent  methods  of  triggering  and  use  of  several 
oscilloscopes  (typically  5)  for  each  shot. 

To  a  lesser  extent  the  following  difficulties  should  be  noted: 

(c)  The  free-surface  velocities  recorded  during  the  experiment 

( typically  15  m/ s)  are  much  smaller  than  those  encountered  in 
large-amplitude  shock  experiments,  but  larger  than  those  in  typical 
acoustical  measurements  in  liquids.  Consequently,  a  novel  form  of 
displacement  interferometer  was  developed  (described  in  Ref.  39)  to 
facilitate  the  velocity  measurements.  There  was  occasional 
difficulty  in  aligning  the  interferometer  and  adjusting  the 
oscilloscopes  because  of  ambient  vibrations. 

(d)  At  the  time  the  experiments  were  performed  (1982),  the  gun  facility 
used  was  a  shared  facility) .  Consequently,  each  shot  required 
special  adjustment  of  the  gun  and/ or  oscilloscopes  for  this 
particular  experiment. 

Technical  suggestions  concerning  possible  future  experiments  of  this  type  are 
noted  in  the  Ad-Hoc  Report.  We  anticipate  issuing  and/or  publishing  a 
complete  description  of  this  research. 
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(abstract) 

Far-field  scattering  patterns  were  photographed  for  oblate  spheroidal 
water  drops  illuminated  by  a  laser  beam.  Near  the  angular  region  of  the 
primary  rainbow,  the  patterns  manifest  the  angular  dependence  of  a  hyperbolic- 
umbilic  (D^*)  diffraction  catastrophe  where  the  control  parameter,  commonly 
denoted  as  C^,  varies  with  the  drop's  shape.  Focusing  corresponding  to 
*  0  was  observed  when  the  drop's  axis  ratio  =  1.3.  The  scattering  is 
novel  both  as  rainbow  phenomenology  and  as  a  manifestation  of  a  D* 
diffraction  catastrophe.  An  ultrasonic  wave  was  used  both  to  levitate  the  mm 
sized  drops  and  to  control  their  shape. 


The  characterization  of  short-wavelength  scattering  phenomena  has  been 

1  2 

advanced  by  the  study  of  diffraction  catastrophes.  *  Examples  include 

1-5  g 

optical  phenomena  and  molecular  collisions.  The  most  familiar  example  of 

7-9 

a  diffraction  catastrophe  is  the  ordinary  rainbow.  The  angular  variation 

of  the  scattering  from  a  spherical  drop  locally  corresponds  to  that  of  a  fold 

1  2 

diffraction  catastrophe,  5  a  result  most  clearly  seen  with  monochromatic 

Q  Q 

illumination.  Uniform  approximations  to  rainbow  scattering  ’  have  also  been 

developed;  however,  these  are  not  essential  to  the  present  research. 

We  studied  the  scattering  of  drops  whose  shape  closely  approximated 

that  of  an  oblate  spheroid  with  the  short  axis,  the  symmetry  axis,  being 

vertical.  The  drops  were  illuminated  by  a  horizontally  propagating  gaussian 

beam  with  a  wavelength  X  =*■  633  nm.  These  drops  were  observed  to  scatter  in 

the  horizontal  rainbow  region  with  patterns  like  those  of  hyperbolic-umbilic 

(classification  D?)  diffraction  catastrophes  . 1* 5 ’ 6  Visible  D+  diffraction 
4  4 

patterns  observed  previously  include  light  transmitted  by  frosted  glass 
surfaces1*"*  and  by  liquid  lenses  clinging  to  tilted  glass  plates. lj 4  Our 
observations  appear  to  be  the  first  report  of  an  umbilic  catastrophe 
associated  with  traditional  rainbow  scattering .  They  suggest  novel  phenomena 
to  be  observed  in  natural  rainbows  or  in  the  evaluation  of  the  partial-wave 
series1'*'  for  scattering  from  spheroids. 

The  ray  diagram  of  the  primary  rainbow  of  a  spherical  drop  is  shown  in 
Fig.  1.  For  the  once-reflected  (twice-refracted)  rays  there  is  a  minimum 
scattering  angle  of  0R  -  138°  and  the  associated  ray  is  known  as  the 
Descartes  ray.  When  the  scattering  angle  9  is  between  0„  and  166°,  there 
are  two  rays  of  the  class  shown.  The  rays  labeled  1  and  2  are  for  0  =  152°. 
In  the  two— ray  region,  intensity  oscillations  produce  the  supernumerary 

Figure  1  is  also  applicable  to  rays  which  lie  in  the  horizontal 


arcs . 
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Descartes 

Ray 


Fig.  1.  Rays  through  either  a  spherical  drop  or  through  a  spheroidal  drop  in 
Che  equatorial  plane. 

Camera 
(Focused  on 
the  Drop) 

Laser  Beam  "  S 

- 7^ 

Camera 
(Focused  at 
Infinity) 

Diffuse  White  Light 

Fig.  2.  Schematic  view  from  above  of  the  apparatus.  The  drop  is  levitated 
11  mm  above  the  plate  and  the  optical  axes  of  the  cameras  lie  near'  a 

horizontal  plane  containing  the  drop.  The  diffuse  light  is  usually  turned 
off. 
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Fig.  3.  Photographs  of  an  acoustical^  levitated  drop  having  a  diameter 
D  -  2.3  mm.  The  drop  was  illuminated  both  with  diffuse  light  and  the  laser 
beam.  The  4  bright  patches  on  the  left  are  due  to  refracted  laser  light. 


Zero-Ray 
Region 


Ca) 


Two-Ray 

Region 


Four-Ray 

Region 


(b) 


T 

H 

1 


Rays  1 , 2, 3, 4 

Fig  •  5.  (a)  Partitioning  of  the  angular  regions  implied  by  the  theory  of 

hyper bolic-umbilic  diffraction  catastrophes  for  sections  having  t  0.  (b) 

Ray  locations  evident  from  the  bright  patches  in  Fig.  3. 
aperture  was  in  the  four-ray  region. 


The  camera's 


equatorial  plane  P  of  the  spheroidal  drops;  ?  cuts  the  dropfs  surface  in  a 

circle  and  both  the  incident  ray  and  the  local  surface  normals  lie  in  P. 

Figure  2  diagrams  our  apparatus  *  A  horizontal  metal  plate  was  coupled 

to  a  piezoelectric  vibrator  which  was  driven  at  a  frequency  of  27  kHz.  A 

curved  reflector  (not  shown)  was  placed  above  the  plate1 s  center  so  as  to 

establish  an  acoustic  standing  wave  in  the  air  below  the  reflector .  Drops  of 

distilled  water  could  be  levitated  in  this  wave  by  an  upward  directed  acoustic 

radiation  pressure.  The  levitator  is  similar  in  principle  to  one  described  in 

Ref.  12.  Figure  3  shows  a  levitated  drop.  The  oblate  shape  is  due  to  the 

13 

spatial  distribution  of  the  radiation  pressure.  In  addition  to  the 

deformation  evident  in  Fig.  3,  the  drop  vibrates  at  27  kHz  with  surface 

displacements  which  are  too  small  to  affect  the  light  scattering.  The  dropTs 

diameter  will  be  denoted  by  D  in  the  plane  P,  and  by  H  along  the  vertical 

symmetry  axis.  Raising  the  ultrasonic  amplitude  increases  D/H. 

The  scattering  was  photographed  for  randomly  polarized  illumination 

from  a  He-Ne  laser.  The  camera  was  focused  on  infinity  such  that  the  pattern 

recorded  was  equivalent  to  that  in  the  far  field  (i.e.  for  distances 

2 

r  »  D  A  from  the  drop).  Horizontal  and  vertical  coordinates  in  each 
photograph  were  linear  in  the  horizontal  and  vertical  scattering  angles  0 
and  3,  where  3  is  relative  to  the  plane  P.  The  camera  viewed  the  region: 
136°  <  0  <  148°,  -6°  <  g  <  6°. 

Figure  4(a)  shows  the  scattering  from  a  slightly  nonspherical  drop 

having  D  s  1.9  mm.  The  intensity  modulations  are  similar  to  other  photographs 

14 

of  supernumerary  arcs  with  laser  illumination.  Within  experimental 

7  9 

uncertainties,  their  spacing  is  as  predicted  by  Airy  theory.  *  Figures  4(b)- 
(f)  are  the  patterns  for  a  different,  more  nonspherical,  drop.  For  Fig.  4(b) 
the  drop  had  D  =  1.39  mm  and  D/H  -  1.23;  somewhat  smaller  values  of  D/H  gave 
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Fig.  4.  Photographs  of  rainbow  region  scattering  patterns  arranged  in  order 
of  increasing  drop  axis  ratio  D/H.  The  horizontal  scattering  angle  0 
increases  from  left  to  right.  The  cusp  in  (b)  moved  into  view  from  the  right 
as  D/H  was  increased,  (d)  is  the  focal  section  of  a  hyperbolic-umbilic 
diffraction  catastrophe  while  (c)  and  (e)  manifest  the  unfolding  of  the 
catastrophe  in  response  to  small  changes  in  D/H. 


patterns  similar  to  Fig.  4(a).  Changes  in  the  pattern  with  changes  in  D/H 

were  also  viewed  through  the  camera’s  lens  by  eye. 

Figure  4(d)  manifests  the  parallelogram  and  caustic  structures 

characteristic  of  focal  section  patterns  of  diffraction  catastrophes. 

1  2 

Consequently,  the  normal  form  9  of  the  associated  diffraction  integral  has 

the  control  parameter  C^  ■  0.  Evidently  the  transverse  control  parameters 

C1  and  C-  are  locally  related  to  g  and  0  -  0  by  a  linear 

transformation.  The  parallelogram  structure  is  sheared  with  an  apex  angle 

-  43.5  +  1°.  Figures  4(b),  (c),  (e),  and  (f)  correspond  to  sheared  and  scaled 
sections  of  D*  diffraction 

patterns  having  C^  /  0.  Inspection  of  the  patterns  of  the  elementary 

14  . 

catastrophes  *  confirms  that  only  the  describes  our  observations. 

When  C^  ^  0,  the  pattern  may  be  partitioned  into  three  angular 

regions  according  to  the  number  of  stationary  points  of  the  phase  function 

within  the  diffraction  integral ^  The  number  of  points  corresponds  to  the 

number  of  once-reflected  (twice  refracted)  rays  which  contribute  to  the 

scattering.  This  partitioning  is  illustrated  in  Fig.  5(a)  for  the  pattern 

shown  in  Fig.  4(c).  That  the  right-most  region  has  four  rays  was  confirmed  by 

direct  observation.  A  camera  having  a  small  aperture,  positioned  as  shown  in 

Fig.  2,  was  focused  on  the  drop.  Figure  5(b)  is  drawn  from  a  representative 

image  so  as  to  show  the  locations  of  the  four  observed  rays.  Rays  1  and  2  lie 

in  the  plane  P  and  correspond  to  those  in  Fig.  1.  Rays  3  and  4  are 

horizontally  directed,  having  g  ^  0°  when  outside  the  drop.  They  are  skew 

rays  since  inside  the  drop  they  do  not  lie  in  a  horizontal  plane.  These  rays 

14  15 

were  not  anticipated  from  previous  studies  *  of  rainbows  of  spheroidal 
drops  which  considered  only  rays  confined  to  a  plane. 

Other  features  observed  were  consistent  with  those  expected  of  a 
diffraction  catastrophe  at  short  wavelengths  X  «  D.  The  axis  ratio  D/H 


which  gives  *  0  should  be  independent  of  D  for  spheroidal  drops  having 
the  same  refractive  index.  To  check  this  prediction,  drops  were  illuminated 
with  diffuse  white  light  (see  Fig.  2)  and  their  profiles  were  photographed. 
Prior  to  each  photograph,  D/H  was  adjusted  so  that  the  laser  scattering 
pattern  was  a  focal  section  as  in  Fig.  4(d).  Photographs  of  six  drops  having 
D  of  1*5  to  2.4  mm  gave  axis  ratios  not  correlated  with  D.  The  average 
measured  D/H  of  1.305  +  0.016  is  a  lower  limit  of  the  true  axis  ratio  since 
any  misalignment  of  the  camera1 s  optical  axis  from  the  equatorial  plane  would 
give  a  systematic  reduction  in  the  measured  D/H. 

1  7 

For  scattering  from  large  spherical  drops,  Airy's  theory  ’  predicts 

that  when  9  =»  8  ,  the  intensity  I  at  a  distance  r  has 
K 

I  «  (D/A)  ^(D /r)^.  Examination  of  the  D*"  diffraction  integral,1 9  ^ 5  6  gives 
2/3  2 

I  «  (D/A)  (D/r)  when  and  9  — all  vanish.  Visual  inspections 

and  photographs  confirm  that  the  scattering  near  the  apex  of  the  focal  section 
appears  brighter  than  that  near  0^  of  the  other  sections  such  as  Fig.  4(b) 
and  (c).  This  focusing  for  ~  0  may  be  observable  in  nature  when  the  sun 
is  close  to  the  horizon.  Raindrops  having  D  -  4  mm  have  an  axis  ratio  -  1.3. 
However,  such  large  drops  are  flattened  below  the  waistline16  and  the  effect 
on  the  focusing  is  not  known. 

These  experiments  were  carried  out  while  P.  Marston  was  at  Jet 
Propulsion  Laboratory  on  a  sabbatical  leave.  The  research  was  partially 
supported  by  the  U.S.  National  Aeronautics  and  Space  Administration  and  by  the 
Office  of  Naval  Research. 
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